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' The anomalous Hall eflFect in disordered band ferromagnets is considered in the framework of 

' quantum transport theory. A microscopic model of electrons in a random potential of identical 

. impurities including spin-orbit coupling is used. The Hall conductivity is calculated from the Kubo 

' formula for both, the skew scattering and the side-jump mechanisms. The recently discussed Berry 

' phase induced Hall current is also evaluated within the model. The effect of strong impurity scatter- 
, ing is analyzed and it is found to affect the ratio of the non-diagonal (Hall) and diagonal components 

■ of the conductivity as well as the relative importance of different mechanisms. 

T— I ■ 

I r I. INTRODUCTION 

a : 

I ' The Hall effect is known to occur in conductors subject to a magnetic field. On a classical level it is explained by the 
^ Lorentz force acting on the charge carriers in a magnetic field. However, an external (or internal) magnetic field is 
not necessary for a Hall response to exist. It follows from basic principles of statistical mechanics that the Hall effect 

1 t ^ can appear whenever time reversal invariance is broken. Indeed it was recognized in the 1950's that a Hall effect 

■ should exist in ferromagnetic metals even in the absence of an external magnetic field, if the magnetic polarization 
of the spin system is coupled to the orbital motion by spin-orbit coupling. An anomalous Hall effect in ferromagnets 

I I has been observed in many systems (for example see^'^). The direct coupling of the moving charge carriers to the 
^ ' magnetic field generated by the spins is much too small to explain the experimental observations. 

O I In a pioneering paper, Karplus and Luttinger'^ worked out a theory of this effect, in which they pointed out the 
existence of an additional term in the velocity operator proportional to the gradient of any electrical potential acting 
on the carriers and to a term acting like a magnetic field (not the dipolar field) . The latter has been identified recently 
as a Berry phase term^"^. It generates a Hall current in equilibrium, and should exist even in the absence of any 
impurity scattering. The Hall conductivity derived from this mechanism is proportional to the Berry phase curvature 
averaged over all occupied conduction band states. The precise dependence of this quantity on the ferromagnetic 
polarization and on the spin-orbit interaction depends on the details of the system considered. 
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' At about the same time, in 1955, Smit^ described a different mechanism known as skew scattering. It is based on 
the fact that electrons in a plane perpendicular to the magnetic polarization scatter from an impurity potential in an 
asymmetric fashion, if they feel the polarization via the spin-orbit interaction. This mechanism has been worked out 
in great detail. It yields a Hall conductivity approximately proportional to the longitudinal conductivity (which is 
1 governed by impurity scattering at low temperatures), to the ferromagnetic polarization and to the spin-orbit coupling 
g ' (assumed to be weak). The relation of this to the earlier works was discussed in^. 

Yet another mechanism for the anomalous Hall effect was proposed by Berger^", the "side-jump" mechanism. It is 
based on the observation that the trajectory of an electron scattering off an impurity is shifted sidewise by the action 
Q ' of the spin-orbit coupling in the presence of a ferromagnetic (or antiferromagnetic) spin polarization. The effect gives 
^ . rise to a Hall conductivity independent of the density of impurities, i.e. of the mean free path. The characteristic 
^ ' length replacing the mean free path is the shift of the trajectory, which may be estimated to be of the order of the 
. ^ , lattice spacing. Hence this contribution is small compared to the skew scattering term , except in the case of short 
mean free path. Although this mechanism is similar to the Berry phase mechanism, and may be shown to originate 
from the extra term in the velocity operator, it involves the nonequilibrium quasiparticle distribution. 

Quantum corrections to the anomalous Hall conductivity have not received much attention so far. The weak local- 
ization correction is cut off by the spin-orbit as well as the phase relaxation rates. Nonetheless, it has been found to 
be of order unity in the disorder parameter (fci?/)~^ , within the skew scattering mechanism^^'^^. In the case of the 
side-jump mechanism, weak localization corrections have been found to be negligibly small^^. Interaction corrections 
have been shown to be absent within the skew scattering model (neglecting Hartree terms) in^^, in the limit of weak 
impurity scattering. For strong impurity scattering, however, a finite interaction correction appears^^, in accordance 
with experimental observation^'^. 

In this paper we will review these different mechanisms of the anomalous Hall effect from a common perspective, such 
that their dependences on parameters are displayed and their relative magnitudes are estimated. We will limit our 
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discussion to two-dimensional or quasi two-dimensional disordered metallic band ferromagnets. As a convenient and 
not unrealistic model of disorder we assume identical short range impurity potentials at random positions, including 
spin-orbit interaction induced by the impurity potential. The strengths of the impurity potential and of the spin-orbit 
coupling will be left as free parameters. In particular we will be interested in strong impurity scattering, which has 
not been considered previously in this context within such a model, to our knowledge. 

II. THE MODEL 

We consider ferromagnetic metallic films with conduction electrons occupying a spin-split band. Transport at low 
temperatures is governed by impurity scattering. We will model the disorder by assuming identical impurities of 
density riimp , at random positions . Electron-electron interaction effects will be neglected. Spin-orbit interaction 
at the impurities will give rise to an anomalous Hall effect, as pointed out in the early papers by Smit^ and Luttinger^. 
This so-called skew scattering arises because the finite magnetization M of the conduction electrons introduces a 
sense of rotation about the direction of M . In addition there is a "side-jump " contribution {^^) , caused by a sideways 
shift of the scattering wave packet due to the spin orbit interaction. The spin-orbit interaction is a relativistic effect 
and therefore rather small for transition metal atoms. As pointed out in^", the mixing of different d-orbitals provides 
a renormalization of the coupling constant leading to an enhancement by a factor of 10^ . We will take this effect 
into account by employing a phenomenological coupling constant Qa of order unity. The single particle Hamiltonian 
of a conduction electron in a ferromagnetic disordered metal , including spin-orbit interaction induced by the disorder 
potential Vdisi^), is given in its simplest form by 



Hi = [-— + Vdis{r)]5,,, - M,T^^, - i{g./A7rn,)[T„„, ■ {Wais x V)], (1) 

where Ua- = {k'^^/^T^) and kpa are the density of conduction electrons and the Fermi wave vector of spin (J , T(7(T/ is 

the vector of Pauli matrices, and g„ is a dimensionless spin-orbit coupling constant. The bare coupling constant is 
given by (/ct"' = {A-Knc){2mc)~^ with m the electron mass and c the velocity of light . Note that c/^ is in general spin 
dependent. We use units such that is half the Zeeman energy splitting of the conduction electron energies caused 
by the ferromagnetic polarization M. Here M is assumed to be oriented along z, perpendicular to the layer . The 
disordered potential will be modelled as Vdis(r) = Sj y(r — R^). We will later average over the impurity positions 
Rj. 

The matrix elements of in the plane wave representation are given by 

{\lW\Hi |ka) = (^— - Maj (5kk/<5aa/ + ^('^ ~ exp[i(k - k/) • Rj]{(5aa/ - igaTaa/ ■ (kxk/)} (2) 

j 

where y(k — k/) is the Fourier transform of the single impurity potential, and k = k/|k| . 

The many-body Hamiltonian is given in terms of electron creation and annihilation operators cj^, Cko- for Bloch states 
I kcr > as 

ko- 

+ Xy(k-k/)exp[i(k-k/)- Rj]{5<,<,, - i5<,w (kxk/)}c+^,Cka, (3) 

k(T,k/(T/ 3 
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III. NORMAL AND ANOMALOUS CONDUCTIVITY IN THE LIMIT OF WEAK IMPURITY 

SCATTERING 



A. Kubo formula and single particle Green's function 

The conductivity will be calculated from the current response functions Lap (^m) by employing the Kubo formula 

1 

= e2 lim — i„;3(iO„) , L^ffiinm) = / rfre'"'"^(r^[ia(T)j>(0)]) (4) 

S2m— >OlZm Jo 

where ja (t) are the components of the current density in the Heisenberg representation , = 2nTm are bosonic 
Matsubara frequencies and the angular brackets denote the thermal average. We will calculate the current correlation 
function within diagrammatic perturbation theory in the impurity potential, averaging over the random impurity 
positions. To simplify the calculations, we will drop the momentum dependence of the single impurity potential, 
y(k) =Vo . The momentum dependence of the spin-orbit term in the potential will be kept,of course, as it is the 
source of the anomalous Hall effect. 

The single particle Green's function GkaiioJn) is defined in terms of the self-energy SkCT(ia;„) as 

Here w„ = nT{2n + 1) are fermionic Matsubara frequencies, and Eko- = £k — -M^cr are the Bloch energies. The Bloch 
states are filled up to the Fermi energy epa for each subband. 

In lowest order in the impurity potential, assuming all momenta to lie in the x-y-plane and after averaging over the 
random positions of the impurities we have 

Sk(7(«W„) = tlimpVo ^{1 - ig„T^„ ■ (kxk/)_^}^ Gklaii^^n) (6) 
k/ 

Since the Green's function as a function of £kcr is strongly peaked at the Fermi energy, one may separate the k/- 
summation into an integral over energy, taking the momenta, in particular that of T,i^cr{i^n) to be at the Fermi 
surface, and an integral over the angle (p/ formed by k/ and the x-axis (note k/ = (cos ip/, sin p/). As a result, one finds 

Skcr(iw„) = -isign{LOn)nimp{'!rNc,)~^[wcr + 2ua-] = -isign{u)n){2Ta)~^ (7) 

where N„ is the conduction electron density of states at the Fermi level for spin subband a and To- are the single 
particle relaxation times. We have defined the dimensionless coupling constants for potential scattering and spin-orbit 
scattering 

= {irN^Vof and = (ga/'^f ■ (8) 



B. Longitudinal conductivity 



The diagonal elements of the conductivity tensor are given in lowest order within our model by the bubble diagram. 
Assuming isotropic scattering potential, there are no vertex corrections in lowest order. We have 

Laa{i^m) = 7" ^ ^ t^k(7,a<^ka (iw„)Gka («W„ - iflrn) (9) 
k,cr 

where Uk<7,a = deka/dka = ka/nicr is the particle velocity, where ma is the effective mass of quasiparticles with spin 

(7 (wc assume parabolic bands) . The integration over energy Sk is only finite and equal to 2'kTcNc if the two poles 
of the G's are on opposite side of the real axis, yielding the restriction on the frequency summation < a;„ < Q.m i 



(a) (b) 

FIG. 1: lowest order diagrams for (a) skew scattering and (b) side-jump contributions. As described in the text, there are two 
diagrams of type (a) and four of type (b). 



and hence (Om/27rT) (identical) terms. The angular average over the velocity factors yields (w^^ ^) = ^Vp^ . As a 
result one finds the sum of Drude conductivities for both spin orientations 

=e' E ^ =e2E^^"'^'^ = e'[4---pKf]-'E-l^./(l + k)' (10) 

filer ^ 
cr a a 

where -0^°'' = ^Vp^Ta- is the diffusion constant. Thus, for weak scattering the conductivity is found to be inversely 
proportional to (1) the density of impurities, (2) the potential scattering cross section Vq , and (3) the factor (1 + ^g'^) 
. Note that D^^Na = efctTct/Stt 1 , considering that in the parabolic band approximation Na = mo-/27r . 



C. Anomalous Hall efFect: Skew scattering contribution 



The lowest order diagram contributing to the Hall conductivity is a bubble with three scattering processes at the 
same impurity, denoted by lines running across (vertex correction needed to give finite angular averages), shown in 
Fig. la, 

l^n k,k/,(T 

Vkk//Gk//CT(iw„)Fk/,k/Gk/CT(ia;„)Gk/(T(it^n - «^^m)14:/k (11) 

k// 



and a contribution with upper and lower line interchanged. Here we use Vk(T,k/(T/ = Vq{1 — ig^T^^ ■ (kxk/)^}. 

Again the energy integrations on Sk, £ki provide the restriction < a;„ < and yield a factor {2'!tt„N„)'^ , while 
the energy integration on the single Gkf/cr{iiOn) gives —isign{u)n)7rN„ and the two angular integrations are finite only 
for the cross tcirms in the product, using {kak/}) = ^Sa/^. The skew scattering contribution to the Hall conductivity 
is then given by 



where in the last equality we used the definition of and of N^- in the parabolic band approximation. It is 

interesting to note that the Hall conductivity is proportional to the spin-orbit coupling for small , while for large 
spin-orbit coupling, gi^ ;» 1 , it is seen to decrease as ■ This regime may be reached only for sufficiently weak 
scattering, such that the conditions for the validity of the above weak coupling calculation , w^, 4; 1 , are satisfied. 
In the absence of spin splitting , i.e. without spontaneous ferromagnetic polarization , the factor r %^ causes a^y^^' 
to vanish. In the limit of weak polarization M (along the z-axis, perpendicular to the ferromagnetic layer), alty^^ is 

proportional to ■ The sign of <J^y*'°'' is given by the sign of (-M^) , provided D^j \n^ , go- are increasing functions 
of the chemical potential, which would be the case for a parabolic energy spectrum. The ratio of anomalous Hall 
conductivity and longitudinal conductivity is small, proportional to (1) the ferromagnetic Zeeman energy in units 
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of the Fermi energy {Mz/sp) (2) the spin-orbit couphng constant (for smaU couphng) (3) the dimensionless 
impurity potential (VqA^q) i where A^o is the spin averaged density of states at the Fermi energy: 

<rV4°2 ^ (^) ^oATo {go). (13) 



D. Anomalous Hall effect: Side-jump contribution 



This effect may be calculated in a straightforward way^^'^^ by observing that the side-jump leads to an additional 
term in the particle velocity due to the spin-orbit interaction. Indeed, the quantum mechanical velocity obtained from 
the Heisenberg equation of motion for the position operator has two terms, 



V = — r = -i[T,Hi\ = ^ + {gal'^T^n„){Taa^yVdis) . (14) 

The matrix elements of v are given by 

(k/a/| V I ka) = ^5kk,^aa/ - i{ga/^TTna) ^ V(k - k/)e^(''-'")-^^ { w x (k - k/)} (15) 

In lowest order in the impurity scattering, the side-jump contribution to the Hall conductivity a^y is calculated 
from four diagrams with an impurity correlator line ended at one of the current vertices and the upper or lower line, 
respectively (see Fig. lb). Since they give identical contributions, we have to consider only the first one. In 2d systems 
confined to the x-y plane, and assuming the magnetization oriented in the z-direction, and the external frequency 
to be positive, 



^ X] X]*^k/a(w„)Gka(w„)Gka(<^n - f^m) X 
Wn k,k/,cr 

X {-ig,Vo'/SF,)[T„„,X^^],^ (16) 

Performing the energy integrations and the angular integrations as in the above, using the definition of t„ and 
multiplying the above by a factor of four to account for all the diagrams, it follows that the Hall conductivity is 
obtained as 

< = (eV27r)$][W(l + j5^)]r^. (17) 



Remarkably, crj|, is independent of the impurity concentration. Since the effective mean free path characteristic of the 
side-jump contribution is rather short, of order g^/kp , it will be important only for dirty samples, when the skew 
scattering contribution is also small. A change in sign of the skew scattering contribution has been found for high 
impurity concentration, in the framework of the Coherent Potential Approximation^^. 



E. Hall current in the clean limit 



The extra term in the velocity operator derived in Eq. (14) involves any potential V{v) acting on the conduction 
electrons. In section 3.4 the potential considered was that due to impurities, Vnsir). In the pesence of an applied 

electric field E, putting 'WV in Eq. (14) equal to — eE, the extra term in the velocity is directed orthogonal to E and 
to the magnetization M, and hence may give rise to a Hall current. Since the velocity operator in this case is itself 
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proportional to the applied field, to obtain the linear response current it is sufficient to average it with the equilibrium 
statistical operator resulting in 



kcr 



where /ko- is the Fermi function. With the help of /kcr = one finds the contribution to the Hall conductivity 
even in the absence of impurities ("clean limit"): 



<y = i E 9. rl^ ■ (19) 

This current flows in equilibrium and does not require the redistribution of particles from excited states into the 
equilibrium state. It is therefore similar to a ballistie current. 

It is a single particle current and it is therefore not stabilized or proteced by collective effects. For that reason one 
may wonder how strongly this current decays as a consequence of inelastic or dephasing process. This will be subject 
of future investigations. 

Recently a Hall current in a perfect crystal lattice of rather similar form has been discussed'*^^ It arises d\ie to a Berry 
phase acquired by electrons moving in the periodic potential of the crystal. A Bloch electron moving in reciprocal 
space under the influence of the combined effect of spin-orbit interaction and a ferromagnetic polarization along a 
path C acquires a Berry phase 

X(k) = - r dkl- X(k/) . (20) 

J c 

Here the Berry vector potential X(k) (of dimension length, and thus describing a shift of the Wannier coordinate of 
the Bloch states within the unit cell), is given by 

X(k) = / <k(r) iVfe u„k(r) . (21) 

icell 

Associated with X(k) is a "magnetic field" 

n(k) = Vfe X X(k) (22) 

acting in k-space. 

The quasi-classical dynamics of Bloch electrons including the Berry phase may be derived from the Bloch Hamiltonian 

ifk = niVfc+Xfc) + e„<,(k) , (23) 

where e„(k) are the Bloch energies (including the spin-orbit interaction and the ferromagnetic polarization), and V{r) 

is the applied external potential. Putting W{r) = — eE, the quasiclassical equations of motion derived from this 
effective Hamiltonian are 

k = eE + ef X B 

r = Vke„k - eE X n (24) 
The additional term in the velocity r leads to a Hall current 

jff = -e^n (rj) X E ^ a^y = e^n (n,) (25) 

where {SI) = n~^Yli\s.a ^<^(l*)/(€k(T) is the average of the Berry magnetic field over all occupied states in k-space 
and n = S^^n^. This average is zero unless time reversal symemtry is broken, e.g. in a ferromagnet, and the spin 
polarization is coupled to the orbital motion. This Hall current has been computed for ferromagnetic semiconductors^, 
and found to be of the magnitude observed in such systems. 

It would be instructive to compare the two types of Hall currents discussed in this subsection within the same model. 
The Berry phase mechanism is derived from the topological properties of Bloch states in a priodic potential, moving 
in k-space under the influence of an external field. The mechanism discussed first does not require the presence of a 
periodic potential, but derives solely from the extra term in the velocity operator due to spin-orbit interaction. 
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IV. STRONG IMPURITY SCATTERING 



The scattering potential of a single impurity is not necessarily small compared to the Fermi energy. If the product 
N^Vo is of order unity, or even larger than unity, or, in other words, if the coupling constants w^, w<t are not small 
compared to unity, repeated multiple scattering off the same impurity becomes important. In technical terms, the low 
order (Born) transition amplitudes we have used in the above must be replaced by the full scattering amplitude. As 
may be expected, this results in a substantial renormalization of the conductivities. 

A. Single particle scattering amplitude 

The repeated scattering of an electron off a single impurity may be described in terms of the dimensionless scattering 
amplitude 

/kcr,k/cr/ = wNoSoadVkaM"^, + VGV + VGVGV + ....] (26) 



where G is the single particle Green's function. V is the bare interaction with one impurity at R = 0. In the following 
we assume short-ranged interaction and again drop the k-dependence of V{k — k'). The dependence on energy is 
governed by the sharply peaked G's, so that we may neglect the dependence on the magnitude of momentum and 
may put |k| =kF<j , the Fermi momentum for spin orientation a. Then the energy integration in intermediate states 

may be done, 

The remaining integrations on angle may be done using (in 2d and neglecting the angular dependence of y(k — k') 
), (k) = , (£2) = 1^ ((k X ki).(ki X k/),)fe, = -i(k . k/) , ((k • ki)(ki X k/),)k, = i(k X k/)„ where we defined 
(O) = / If O. Within these approximations the potential scattering and spin-orbit scattering terms do not mix. Odd 
terms in V do not depend on s^^ whereas even terms do. The terms of perturbation theory may be summed up to 
infinity as a geometric series, with the result 

II) '2'}i 

/ka,k,a = ^ - (k X k/) — ^ - is^„ [Wo + 2Uo (k • k/)] . (27) 



Here we defined Wc = Wcl{}- -\- w^)-, and Ur, = Ual(\ -\- Ua), where Wo- = (TrA^o-Vo)^ and = {ga/'^Ywo , and 
all quantities depend on the spin orientation a (suppressed in the following, except in the final expressions involving 
spin summation). 

B. Single particle relaxation rate 

The single particle relaxation rate is obtained from the imaginary part of the self energy 

— = 2sc^„ImI]k(7(w„) = 2s^^„niTOpIm/ka,ka = "^-^{wa- + 2Ua) , (28) 



where riimp is the density of impurities (number of impurities per volume). One observes that is proportional to 
the Fermi energy, the average number of impurities per electron and the dimensionless factor {w + 2m) , expressing 
the effective scattering strength per impurity. 



C. Particle-hole propagator 



The particle-hole propagator rkk/(q; f-mf-n — ^m) is an important ingredient of vertex corrections of any kind. Here 
k + q/2, k — q/2 are the initial , k/ + q/2, k/ — q/2 the final momenta and e„, e„ — fim are the Matsubara frequencies 
of the particle and the hole line, respectively. F satisfies the following Bethe-Salpeter equation (we have defined 
dimensionless quantities r,t by multiplying both with a factor {2nNoTo)) 
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rkk/(q; icn, i^n - i^m) = tkk/(q; icn, i^n - i^m) + {^irNaTa) ^ 

^ ^ ^ ^kki (q; i£ni i^n ~ i^m.)Gki+ci/2,a{i^n) 



ki 



(29) 



The kernel of this equation is the impurity averaged particle-hole scattering amplitude (we consider only the case of 
equal spin of particle and hole) 



ikk/(q; o-;k/+q/2,o-(ien)/k/- q/2,(T;k—q/2, (7(^^71 ^^m) 



(30) 



It is useful to represent the operator tkk/(q = 0) in terms of its eigenvalues . Assuming isotropic band structure, 
the eigenfunctions Xm(k) = cxp{irrnp) are those of the angular momentum operator component L^. The operator 
*kk/ (1 ~ ^) represented as 



*^k;(q = 0) = $]A„Xm(k)x™(k/) and t^+i<i = 0) = [t+^,{ci = 0)]* 

m 

The eigenvalues for Sp = sign(e„) = 1, Sh = sign(en — f^m) = —1 are given by 

Ao = l , Xi =2wu{w + 2u)-'^{l+iSp-^T^J 



(31) 



(32) 



A2 = — (w + 2u) ^ {u — 1 + 2iSpy/uT^^) and \-m = Kn- 



(33) 



The energy integral over the product of Green's functions in the integral equation for Fkk/ may be done first, after 
expanding the G's in Q,m and q, where q • Vj. = qvp^S.-^) ■ Expanding Fkk/ and tkk/ in terms of eigenfunctions Xm(k) 
, Tkk/ = Emr„jk/Xm(k), one obtains 



r^k/ = A„jXm(k/) + A„j| l-r(|0„|+Dog^) - ^vpqr Tm-vkiXi{^} + 

+ r„+ik/xi(q) -^{vpqrf r„_2k/X2(q) + rm+2k/X2(q)]} 



(34) 



The case rn = needs special consideration, because particle number conservation causes Fok/ to have a pole in the 
limit Q.m, 5' — » , here expressed by Aq = 1 . 

The complete particle-hole propagator in the regime vpqr < 1 is given by 



" 10 l^i^n^2 ^k7k/ + ^ A,„Xm(k)Xm(k/) ; A„ = 

|i'rn| + J^q ^ ,„ J- ■^■n 



(35) 



7k = 1 - ^Vpqr J2 ^mXm(k)Xm(q') = 1 - «T X] '^™Xm(k)(q • Vfe,Xm(k/))k/ 



m=±l 



m=±l 



(36) 



7k = 1 - ir X A„Xm(k)(q • Vfc,Xm(k/))k/ 



m=±l 



(37) 
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where the renormalized diffusion constant is defined as 

D^ = D^° \}~ ^\^ , where Ai = Re Ai. (38) 
|i - All 

The vertex corrections of the current vertices j'ka and j'ka (for the incoming and outgoing current) are obtained by 

jkcr,a(5) = Vka + (Vk/ark/k)k/ 

ika,a(5) = Vka + (Vk/arkk/)k/ (39) 

Note that jka 7^ ( jka)* , as the eigenvalues Xm are complex valued, in general. 

D. Conductivity tensor: Skew scattering 

Multiple scattering off the same impurity generates momentum dependence of the effective scattering amplitude and 
gives rise to vertex corrections to the current vertex as calculated above. The current-current correlator is now given 

by 

Laf3{ii^m) = T^^Vi^cr,aha,0i<l = 0)Gkcr (iWn)Gkcr (iw„ - i^m) (40) 
lii„ k,cr 

Using the explicit expressions for the components of the current vertex 

IkaAQ = 0) = VFa[{l + A'i)4 - Kky] ; jka,j/(g = 0) = VFa[{l + K)ky + Ai%] (41) 

where the terms involving the diffusion pole drop out (they involve an angular average {q^ Qy) = 0), one finds for the 
conductivity tensor 

where we have defined A^ = RcA"^ and A" = ImA^. We recall that = Aict/(1 — Ai^) and A'^^ = 2wu{iij + 2u)~^, 
X'l^ = X'l^-^T^^. Defining the tensor of diffusion coefficients as 

= \vl.r'J , = ^aa[A7(l + A')] = "^^x (43) 

where r*'' = ro-(l + A') is the momentum relaxation time, we may write 

<% = Y,N,Dlp. (44) 



E. Side-jump mechanism 

In the strong scattering regime the bare impurity scattering potential needs to be replaced by the scattering amplitude. 
In addition the vertex corrections to the current density operator have to be applied. There are two new diagrams 
involving the scattering line denoting the spin-orbit term in the velocity operator, Eq. (14), and ending at one of the 
current vertices, framed by two (instead of only one) scattering amplitude lines (see Fig. 2). 
The result of adding the six diagrams is 

a-i = ^ Y t'- \±}^lz (45) 
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F. Limiting cases of Anomalous Hall conductivity and comparison of contributions 



In the limit iVoK 
is then given by 



>> 1, the results simplify considerably. The sum of the skew scattering and side jump contributions 



'xy 



xy 



E 



1 1 



V Tlim. 



(46) 



It is seen that in this case the skew scattering term dominates even in the limit of large impurity concentration. 
In the case of strong potential scattering, but weak spin-orbit interaction, w » 1, but u « 1, we find 



2 '^^riimp) 2TrWa 



(47) 



By comparison, the clean limit Hall conductivity a^y depends only on the spin-orbit coupling gcr (see Eq. (19)) or in 
the case of the Berry phase contribution, on the Berry magnetic field f2 (see Eq. (25)). The former is not dependent 
on impurity scattering, as well as the latter. One observes that the signs of a^y and a^^y are the same. In the weak 
scattering limit cr|-^ is a factor of two larger than a^y, but with increasing scattering strength it drops to values much 
less than a^y. The sign and magnitude of a^y have not been calculated for the model under consideration here. 

Within the model considered here the skew scattering contribution will dominte all other contributions in the limit 
nimp 0. In a more refined model, however, assuming scattering centers with weak (or no) spin orbit interaction of 
density n„ in addition to the skew scattering centers with density rUmp, the mean free path in the limit riimp — * will 
be limited by the normal scattering processes, and consequently a^y x riimp/n'^ and cr*|, oc riimp/nn will tend to zero 
for n 



■%mp 



0. In this case the clean limit contributions al,, and erf;, will survive. 



V. CONCLUSION 



The anomalous Hall effect is a surprisingly rich phenomenon with many interesting facets. It requires broken time 
reversal symmetry as realized in magnetically ordered states and, as far as the symmetry breaking occurs in spin 
space, it requires sufficiently strong spin-orbit coupling. The quantum nature of electron scattering by any impurity 
potential including spin-orbit coupling leads to a right-left asymmetry of the average scattering probability in third 
order of the potential ("skew scattering") and to an extra contribution to the velocity operator ("side-jump" effect). 
Both contribute to the Hall conductivity in zero (or low) magnetic field. In addition, an extra term in the velocity 
arises directly from the applied electric field, for a uniform system as well as in the periodic potential of the crystal. 
The latter contribution has been shown to be a consequence of a Berry phase associated with the motion of Bloch 
electrons in momentum space. The extra term in the velocity proportional to the modulus of the applied electric field 
and directed perpendicular to it has a finite equilibrium expectation value, yielding a Hall current for the clean system. 
A systematic experimental investigation showing the existence of all these contributions in a controlled fashion does 
not exist yet. Further experimental work on well-characterized systems as a function of magnetization, impurity 
concentration, for different strengths of spin-orbit interaction is necessary. Theoretical models for these systems need 
to be refined, using band structure calculations, to model impurity scattering and spin-orbit interaction in a more 
realistic way. 
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